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Abstract 

An effective Lagrangian approach to describe the dynamics of confinement and symme- 
try breaking in the process of quark-gluon to hadron conversion is proposed. The deconfined 
quark and gluon degrees of freedom of the perturbative QCD vacuum are coupled to color 
neutral condensate fields representing the non-perturbative vacuum with broken scale and 
chiral symmetry. As a first application the evolution of gluons emitted by a fragmenting 
high energy qq pair from the perturbative to the non-perturbative regime with confinement 
is studied. For reasonable parameter choice the solution of the equations of motion leads 
to fiux tube configurations with a string tension i ~ 1 GeV/fm. 
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The transformation of color charged quarks and gluons in high energy QCD processes 
into colorless hadrons is commonly believed to be a local phenomenon. The universality 
of parton fragmentation observed in hard processes as e+e" annihilation, deep inelastic ep 
scattering, Drell Yan, etc., is strongly supported by experimental observations. Thus, a 
consistent description of the local hadronization mechanism must be independent of the 
details of the partons prehistory and should apply also to hadron-hadron, hadron- nucleus, 
or nucleus-nucleus collisions. 

To date most of the theoretical tools to study properties of QCD are inadequate to 
describe the dynamics of the transformation from partonic to hadronic degrees of freedom: 
Perturbative techniques are limited to the deconfined, short distance regime of high energy 
partons [1], QCD sumrules [2] and effective low energy models [3] are restricted to the 
long distance domain of hadrons, and lattice QCD [4] lacks the capability of dynamical 
calculations concerning the quark-gluon to hadron conversion. On the other hand, phe- 
nomenological approaches to parton fragmentation [5], are mostly based on hadronization 
models with adhoc prescriptions to simulate hadron formation from parton decays. 

In this letter I will propose a universal approach to the dynamic transition between 
partons and hadrons based on an effective QCD field theory description. The key element is 
to project out the relevant degrees of freedom for each kinematic regime and to embody them 
in an effective QCD Lagrangian which recovers QCD with its scale and chiral symmetry 
properties at high momentum transfer, but yields at low energies the formation of symmetry 
breaking gluon and quark condensates including excitations that represent the physical 
hadrons. I will first formulate the general field theoretical framework. Subsequently I shall 
demonstrate its applicability to the dynamics of parton-hadron conversion by studying 
exemplarily the evolution of gluons produced by a fragmenting quark-antiquark pair. The 
aim of this letter is to give a lucid presentation of the conceptual framework, details of the 
work reported here will be given elsewhere [6]. 

The effective QCD Lagrangian is represented in the form 
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is the QCD Lagrangian in terms of the quark (ipi) and gluon fields (Aa) and F^" = 
d^A^ + gJahcA'^^A'^^. The eff-ective Lagrangian [7] 
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introduces a scalar gluon condensate field % and a meson field U = f^^ exp (i Y^^=o 
for the nine pseudoscalar meson fields (f)j i^f^ = 93 MeV), with non- vanishing vacuum 
expectation values Xo = (0|x|0) 7^ Uq = c(0|C/ + i7^|0) ^ that explicitly break 
scale, respectively chiral symmetry. In (3), c is a constant, rriq = diag(m„, m<j, rrZs) and h 
is related to the conventional bag constant B hy B = hx^l^- The key ingredient in (1) is 
the connection between the scale and chiral symmetric, short distance regime of colored 
fiuctuations and the world of colorless hadrons with broken symmetries. It is mediated by 
the coupling between the fundamental quark and gluon degress and the eff"ective fields % 
and U through coupling functions g{x) i{x)i 

= ^F,^,^Fr - A9ix)i'^■ (4) 

The couplings are chosen [6] in the spirit of the Friedberg-Lee model [8] as g{x) = ^0^/(^ — 1) 
and ^(x) = 9{x)v^{^v — 4) with v = x/Xo which generate absolute confinement of quarks 
and gluons [9], here intimately connected with the formation of the gluon condensate Xo 
and the qq condensate Uq- The non-linear couplings have the properties that ^ = g = for 
X = at short distances, but in the long range regime ^^1, gi^ooasx^ Xo- 

The specific form of AC = C[x, U] + C[^l), A, x] is motivated by the dual QCD vacuum 
picture [8]: High momentum, short distance quark-gluon fiuctuations (the perturbative 
vacuum) are embedded in a confining background field x (the non-perturbative vacuum), 
in which by definition the low momentum, long range fiuctuations are absorbed. In analogy 
to a thermodynamic system in equilibrium with a heat bath there can be a net fiow of energy 
between the system and the heat bath environment such that the bare energy of the system 
is not conserved. However the free energy of the system, here high momentum quarks and 
gluons, is constant [10]. It corresponds to the conserved physical energy momentum tensor 
9 1^1, with the well known trace anomaly [11]: 

0,, = (1 + 7™) E ^^^9 + ^P^'^P"' ' (5) 
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where 7^ is the anomalous mass dimension and (3 the QCD beta function. This anomaly 
constraint is modeled [7] by the third and fourth term in (3) with the correspondence 
C/o = (0|gg|0) and -Bxt = {0\{(3/2g,)F^,F'^''\0). 

The following remarks concerning the effective QCD Lagrangian (1) are important: 
First, at short distances the exact QCD Lagrangian is recovered, since x = U = and 
^ = = 0, whereas the long distance QCD properties emerge as x/Xo^ U/Uq 1 [7] and 
no colored quanta survive. Second, the problem of double counting degrees of freedom has 
to be carefully inspected. Although it does not arise in one-loop calculations (to which 
I will restrict here), processes with e.g. two-gluon exchange could also be contained in 
the exchange of a color singlet x-quantum. A minimal possibility to avoid this problem 
is a rigid separation of high and low momentum modes, in a way that it is exclusive but 
complementary. Third, £[%, U] embodies the correct QCD scaling properties. It accounts 
for the anomaly (5) without which Poincare invariance would be broken, and the mass 
of the proton would come out wrong, since = {p\0'^\p)- Fourth, there is no need for 
renormalization of £[%, U] since % and U are already interpreted as effective, composite 
fields with loop corrections implicitely included. 

The effective QCD field theory defined by (1) is readily applicable to describe the 
dynamic evolution from perturbative to non-perturbative vacuum in high energy processes. 
In the following I shall consider as an example the fragmentation of a qq jet system with its 
emitted bremsstrahlung gluons and describe the evolution of the system as it converts from 
the parton phase to the hadronic phase. A time-like virtual photon in an e+e~ annihilation 
event with large invariant mass ^ A^ is assumed to produce a qq pair which initiates a 
cascade of sequential gluon emissions. The early stage is characterized by emission of "hot" 
gluons far off mass shell in the perturbative vacuum. Subsequent gluon branchings yield 
"cooler" gluons with successively smaller virtualities, until they are within Ql, where Qq is 
of the order of = ibxo — 1 GeV. At this point condensation sets in, or loosely speaking, 
the "cool" gluons are eaten up by the color neutral scalar condensate field % representing 
the non-perturbative vacuum. Restricting the analysis to the gluonic sector, the equations 
of motions that derive from (1) reduce to: 

dA<x)Fn = -gs<x)fa,cA,,,Fr (6) 
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where k = 1 — ^ and V{x) is the potential evident from (3). The key problem is the first 
equation, since the gluon fields drive the dynamics of the x-field which in turn feeds back 
via K,[x)- Note that the [/-field does not couple directly to the gluon field. The procedure 
in the following is therefore to 'solve' eq. (6) as a function of k and then to insert the 
solution into (7), so that one is left (aside from the simple third equation) with a single 
non-linear equation for %. As a simplification I shall work in the mean field approximation, 
i.e. separating off" the quantum ffuctuations of the x-field, % = % + %, where % is a c-number 
and X a quantum operator (similarly for U). 

Solving the equation of motion (6) for the gluon field is equivalent to the calculation 
of the complete Greens function with an arbitrary number of gluons, a formidable task 
indeed. Instead, I will restrict to evaluate the 2-point Greens fuction only (Fig. 1), i.e. 
the full gluon propagator which includes both the gluon self-interaction through real and 
virtual emission and absorption, and the effective interaction with the confining background 
field X [6]- In the framework of "jet calculus" [12], this gluon propagator describes how 
a gluon, produced with an invariant mass Q^, evolves in the variable x and the virtuality 

(or transverse momentum k\) through these interactions. In one-loop approximation 
it is obtained by calculating the corresponding single loop cut diagrams. In the present 
case, one has in addition to the usual gluon branching and fusion processes, g gg and 
gg g, contributions from energy transfer and two gluon annihilation processes g gx 
and gg %, respectively. This is illustrated in Fig. 1. 

The determining equation for the propagator corresponds an evolution equation for the 
gluon distribution, which is defined [13] as average number of gluons at 'light cone time' 
r"*" = in the multi-gluon state \P) with light cone fractions x = k'^ j P'^ in a range dx and 
transverse momenta in a range d^fc^ with respect to the jet axis: 



Performing a standard Mellin transformation to the moment representation g(u),k\) : = 
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dxx"g{x, k\), the evolution equation reads: 
where 7(0;, fc^.) plays the role of a generalized anomalous dimension (A^ = ^(x)^/4x), 
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In Leading Log approximation one finds: 
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Here = N^, p = 0(1) is a parameter, S(uj) = il)(uj + 1) — V'(l) with the Psi (Digamma) 
function '>p{z) = c?[lnr(z)]/c?z and — V'(l) = 1e = 0.5772 the Euler constant. Fig. 2 shows 
7(0;, k\) for diff"erent values of k\. 

The formal solution of eq. (10) is obtained by complex integration in the a;-plane, 
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which yields the a;-distribution of gluons as a function of the evolution variable k\. For 
the full anomalous dimension (11) this inversion must be done numerically [6]. However, 
the main features of the evolution of the gluon distribution in the presence of the x-field 
become already evident when making some simple analytic estimates. For simplicity I will 
divide the kinematic range of the ^^.'^volution into two distinct domains: (i) ~> k\> 
and fci/A^ > 1 so that ~ and 9A^/91nfci ~ 0, (ii) Ql > kl > and fc^/A^ ~ 1 
with A^ (4x)~^ and dX-^j d\a.k\ > 0. In this case, the gluon distribution g(x, fc^.) can be 
evaluated analytically [6] which yields for k\ >, Ql the well known Leading Log solution of 
QCD [14], but for k\ < Ql a substantial suppression at small x is found to set in, refiecting 
the "condensation" of gluons into the confining background field %. The corresponding k\- 
dependence of the total gluon multiplicity Ng{k\) = dx g(x, k\) can also be estimated 



in the above approximation with the result: 
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In the region (i) Q'^ > k\ > Qq, the multiplicity coincides with the QCD result [14], char- 
acterized by a logarithmic growth as the gap between the hard scale and k\ increases. 
On the other hand, in the region (ii) Ql > k\ > A^, the exponent is always negative so that 
the number of gluons rapidly decreases and vanishes at k\ = 0, ensuring that no gluons 
and therefore no color fluctuations exist at distances R >, A~^. This behaviour is evident in 
Fig. 3. 

Returning to eqs. (6)-(8), the coupled system of equations can now be solved numer- 
ically by utilizing the deflnition (9) together with the solution for the gluon distribution 
(10) and (13). Since % is taken to be equal the mean fleld %, the multiplication of eq. (7) 
by (P| and |P) from the left and right, respectively, afl"ects only the Ff^^F'^" term which 
gives the gluon distribution by virtue of (9). With the solution of the gluon distribution 
inserted, the remaining task is to solve eq. (7) for the x-fleld. An interesting phenomeno- 
logical application [15] is to calculate the string constant, which characterizes the linearly 
rising potential between the gg-pair due to the gluon interactions. Considering the frag- 
mentation of a heavy qq pair, an adiabatic (Born-Oppenheimer) treatment is a reasonable 
approximation, because one can view the i*'^ gluon as being emitted from the qq pair plus 
gluons gi,g2i ■ ■ •jfl'i-i, with the spatial coordinates of these "sources" being frozen during 
the emission of the gluon i [16]. Using the separation Rgg of the receding qq as a measure of 
the typical gluon transverse momenta k\ ~ R~? and minimizing the energy per unit length 
of this system, one obtains in the adiabatic approximation for each gluon conflguration at a 
given Rgg the form of % and the string tension t. The total energy per unit length t receives 
contributions from both the x-fleld in the low energy regime and the gluon fleld in the high 
energy domain, which on account of (5) is given by [6], 



'r 



99 



/ 



dA 



v{x) 



4a. 



F,^,^Fr\P), (16) 



where A is the cross-sectional area of the flux tube of the gluons between the qq, perpen- 



dicular to Rgg. In one loop order (3 = — with b = 33/(12x) for Nf = 0. Using eq. (9), 
assuming axial symmetry, and minimizing t with respect to % yields following non-linear 
integro-differential equation (r is the radial coordinate perpendicular to Rqq): 



This equation is to be solved subject to the constraint [15] that the total color electric 
flux through a plane between the q and q equals the color charge on one of them, and the 
boundary conditions x'(0) = 0, x(oo) = Xo- The solutions for t, % and k are shown in Fig. 4 
for the reasonable parameter choice Xo = and bag constant B = (150 MeV)^ [6]. As the 
quarks separate with increasing Rgg, the gluons flrst multiply which results in a growing 
string tension, but then gluon condensation sets in, yielding a saturating behaviour of t 
(Fig. 4a). For Rgg ~ 1 fm the gluon fleld is completely conflned within a flux tube of radius 
r ~ 1 fm (Fig. 4b). 

In perspective, the Lagrangian approach presented here to describe the conflnement 
dynamics of high energy partons in the framework of efl"ective QCD fleld theory has the 
potential to be extended to a systematic quantum description of hadronization. With the 
inclusion of quark degrees and quantum fluctuations of the %- and [/-flelds, one could 
calculate e.g. the mass spectrum of glueball and meson excitations as physical hadrons. 
Ultimately one would like to address the microscopic dynamics in full 6-dimensional phase- 
space, with explicit inclusion of the color flow throughout the evolution. The possible 
applications are manifold. One particular interest is the expected (non-)equilibrium QCD 
phase transition in high energy systems as in heavy ion collisions or the early universe, an 
issue which could be addressed along the lines of Campbell, Ellis and Olive [7]. 
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FIGURE CAPTIONS 



Figure 1: Diagrammatic representation of the two-point Greens function of gluons, in- 
cluding both the gluon (self) interactions and the effective interaction with the confining 
background field % (indicated by the dashed lines). This gluon propagator describes the 
evolution of a gluon from a chosen cascade branch in x and k^, starting from Xq and Q^. 

Figure 2: The anomalous dimension 7(0;, fc^.) of eq. (11) versus u) for diff"erent values of 
kl. 

Figure 3: Evolution of the gluon multiplicity Ng[k]_) from down to A^. At first the 
gluons multiply, but at fc^ ~ a rapid condensation sets which is complete at A^ [Qo = 1 
GeV, A = 0.23 GeV). 

Figure 4: a) String tension t versus separation Rgg of the qq pair, and b) the solutions 
for X and k at Rgg = 1 fm versus r which is the radial coordinate perpendicular to Rgg 
(Xo = Aand 51/4 = 150 MeV). 
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